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SYSTEM TOPOLOGY 


A coupled dynamical system is defined as an assembly of rigid/flexible bodies 
that may be coupled by kinematical connections* Large relative displacement 
and rotation are permitted. The interfaces between bodies are modeled using 
hinges having 0 to 6 degrees of freedom* A hinge is defined as a pair of two 
material points, one on each of two adjoining bodies* A reference body is 
arbitrarily selected and it is assumed for convenience that the reference body 
is connected to an imaginary inertially fixed body. For consistency, a ficti- 
tious hinge is assigned to the reference body by assuming Pq (See Figure), an 
inertial point* Thus the number of hinges equals the number of bodies in the 
system (as shown in Figure 2). 



cut joint 





Inertial Reference Point 
and Fraae 


Figure 1 An Arbitrary Topology 


cut joint 



Figure 2 Equivalent Tree Topology 



Figure 3 A Typical Hinge 



CONSTRAINED DYNAMICAL SYSTEM 


For a mechanical system of n flexible bodies in a topological tree configura- 
tion, the equations of motion are presented in Reference 1. Lagrange's form 
o D Alembert s principle was employed to derive the equations. A detailed 
discussion of the approach is available in References 1 and 2. Equations 
( 4), (15) in Reference 1 are the motion equations for the system of Figure 2. 
These equations are augmented by the kinematical constraint equations. This 
augmentation is accomplished via the method of singular value decomposition 
(see Reference 3). 
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CONSTRAINED DYNAMICAL SYSTEMS 


• LET q 1 .....,q n DEFINE THE CONFIGURATION OF 
UNCONSTRAINED SYSTEM 

• CONSTRAINT EQUATIONS 

A q = B 

• NUMBER OF INDEPENDENT COORDINATES 

= n- RANK A 

• EQUATIONS OF MOTION 

M q = £ + f 

c 

f - FORCES/MOMENTS OF CONSTRAINT 
c 

• PROBLEMS OF PRACTICAL INTEREST 

- SIMPLE NON HOLONOMIC OR HOLONOMIC CONSTRAINTS 


n+ m 

2nd 

ORDER 

D. E. 

OR 2 n + m 

1st 

ORDER 

D. E. 


A 

■ XI 

m < n 
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SINGULAR VALUE DECOMPOSITION 


L mxn MATRIX 

THERE EXIST ORTHOGONAL MATRICES U 

r X I 0 


in 


AND 


SUCH THAT U LV 


r X I 0 1 

— f — 

. 0 I 0 . 


X - D I AG (X 1 >X 2 > >X r ) 


r = RANK L 


Aj - NONZERO EIGENVALUES OF L T L 
L - USV T 


r x I o I r v; i 
' 1 1 1 1 [ ; 


V 2 SPANS NULL SPACE OF L 
U 2 SPANS NULL SPACE OF L 


LV, 


+ r x 1 1 0 1 

• - V — u 
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PSEUDO INVERSE OF L 
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nxn 
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APPLICATION OF SVD TO CONSTRAINED SYSTEMS 


• Aq ■ B 

q * A + B + V 2 i Z COLUMN VECTOR (n-r) 

V 2 Z SOLUTION OF Aq - 0 
A + B PARTICULAR SOLUTION 

• DIFFERENTIATE 
Aq - ft - Aq 
Aq - B* 

q * A + B + V 2 2T 

• Z ' s ARE REDUCED SET OF n-r COORDINATES 
V 2 IS THE DESIRED TRANSFORMATION 
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• Mq = F + F* 

PROJECTION ON NULL SPACE OF A 

VjMq = VjF + VjF* 

SUBSTITUTE FOR q 
V 2** V 2 4 = V 2 F + V 2 F ° “ V 2 MA+B 

• F C - A T X A - LAGRANGE MULTIPLIERS 

v!f c « vJa t a 
2 2 

- [ av 2 ] T X - 0 

• iVjMVjli* = VjF - V*MA + B GOVERNING DIFF. 

q = A + B + V 2 4 


EQN. 


1204 


DAMPING LAV 


• ACTUAL DAMPING LAV: 

F f < t l) “ £ ( LOAD OVER CONTACT AREA,* 

x < t ) . 0 < t < tj. tj, ...) 

F f = FRICTION FORCE 

x(t). 0 < t < tj, REPRESENTS THE HISTORY OF MOTION 

PRIOR TO 

• IN GENERAL THE DAMPING LAV IS VERY COMPLICATED 




COULOMB DAMPING 

F f ( tj) = f (N, S I G N ( x ( tj) ) , 

N - NORMAL LOAD 

|i g - STATIC COEFFICIENT OF 

|i d - DYNAMIC COEFFICIENT OF 

( 1 ittj) > 
0 x(t 1 ) = 

-1 xCtj) < 


S' 


FRICTION 

FRICTION 

0 

0 

0 
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THE MOST COMMON DRY FRICTION DAMPING LAW 

", ' "d 0R F . - F d 

MANY TIME DOMAIN AND FREQUENCEY DOMAIN STUDIES HAVE 
BEEN PERFORMED FOR HARMONICALLY EXCITED SYSTEMS WITH 
DRY FRICTION DAMPING. THERE ARE THREE BASIC METHODS 
OF ANALYSIS 

EXACT 

HARMONIC APPROXIMATIONS 
TIME INTEGRATION 

OUR APPROACH IS TIME INTEGRATION 
STICK/SLIP MOTION 



t Q BECOMES AN UNKNOWN. 


STUCK HINGE RESULTS IN ADDITIONAL CONSTRAINTS ON 
KINEMATICAL VARIABLES. 





COULOMB DAMPER ALGORITHM 



• Y CHANGES SIGN 

COMPARE |Y DIFFERENCEl TO A PRESCRIBED e r*/ 0 

IF GREATER THAN e GO BACK TO PREVIOUS STEP AND REDUCE 
THE STEP SIZE AND REPEAT UNTIL |Y DIFFERENCE! < e 

ACTIVATE THE CONSTRAINT Y = 0 AND COMPUTE CONSTRAINT 
FORCE f c FOR THIS CONSTRAINT 

IF f c OVERCOMES f d SLIP CONDITION. DEACTIVATE Y = 0 

IF NOT, STICK CONDITION, RETAIN Y = 0 AND KEEP 

COMPARING f c WITH f d UNTIL SLIP CONDITION 
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YJI TRANS. AXIS 1 . HINGE 2 c*x # / 

YJI TRANS. AXIS i . HINGE 1 
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TIME (SEC) 




1210 


TIME (SEC) 






YJI TRANS. AXIS 1 . HINGE 1 Oisf « 3 

YJI TRANS. AXIS 1 . HINGE 2 
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TIME (SEC) 
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TIME (SEC) 
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